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PRELIMIiSrARY WORK OIST THE DETERMINATIOISr OF THE 
LAW OF THE PROPAGATION OF HEAT IN THE INTE- 
RIOR OF SOLID BODIES. 

By B. O. Peirce, Jr., and Edward B. Lefavour. 

Presented Oct. 10, 1877. 

For a long time it has seemed probable, as stated in a paper * pub- 
lished last spring in the Proceedings of the Academy, that the flux of 
heat in any direction x, in a solid body, can be written, — 

dAv) 

where v is the temperature, c a constant different for different sub- 
stances, and f(^v) an undetermined function of v. The object of this 
paper is to show that such a function f{y) can be found. 

This is not of necessity possible ; for denoting by X, Y, Zy the com- 
ponents in three rectangular directions of the vector function which 
represents the flux, we have 

dy 

^ dF{x, y, z) _ 



dz 
whence 

— c dF{x, y, 2r) = X 6£r 4- T dy -\- Z dz 

This equation is integrable, and F(x, y, z) can be found only when 
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* Note on the Determination of the Law of Propagation of Heat in Solid 
Bodies. 
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Our experiments were directed to determining whether these condi- 
tions can be satisfied, and 

F{x, y, z) =f{v). 

When a body heated in any way reaches a final state, — that is, a 
state where just the same quantity of heat enters each portion during 
a given time as leaves it, — the function y(y), if it exists, must satisfv 
the equation 

d.YV^) dy'^Av) c/.2/^,;) ^ 

c/t2 "^ c/^2 "T ^^2 

One, and only one, solution of this equation corresponds to each set of 
physical conditions. Since v, and any function of v, are constant along 
the same surfaces, if, when the body is in a final state, v is constant 
along each surface of the family cp (x, y, z) = k, where g) (a?, y, z) is 
the solution of Laplace's equation, corresponding to the given physi- 
cal conditions, then it is always possible to find a function of the tem- 
perature alone, which shall satisfy Laplace's equation as a function of 
X, y, z, or, what is the same thing, shall be equal to cp (x, y, z) through- 
out all space. 

For let u and v be two functions of x, y, z, such that u =z k, and 
V = c represent the same family of surfaces, then denoting by du the 
total differential of u, and by d^u the partial differential relative to x, 
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d»u dyU dgU 

'dF "dii dz . . 

5^ = 5 = 5^ =i/;(a:,2^,^). 

dx dy dz 



If p is any variable. 



And 



dpU dxU dpX 1^ dyU dpy d»u dpZ 

dp dx dp dy dp dz dp 



dPv dxV dpX . dyV dpy d^v dpZ 

dp dx dp "^ dy dp dz dp 
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or, substituting, 

dpii 



dpV 

Similarly the ratio of the correspondiog total differential coefScients is 
du 

Whence, changing the variable and integrating, 
dyii du 

The partial and total differential coefficients of u taken relatively to v 
cannot be equal, if u involve any other quantity than v, 

.•. u =f{v) 

In short, when a body is heated in any manner whatever, there 
must exist a functiony(2;), the same for all bodies, whose derivative in 
any direction, when multiplied by a constant depending on the nature 
of the body, gives the flux of heat in that direction, provided v is found 
constant along the surfaces gp (a?, y, z) == k, which belong to the solu- 
tion of Laplace's equation for that particular case. 

The first case open to direct and satisfactory experiment is where a 
plate of metal is heated at two points, and exposed to the air only at 
its edges. The isothermals in this case belong to the family A log r^ 
+ i? log rg = ^ ; 



This latter form of the equation shows that the two points need not 
be heated equally. The solution for three dimensions cannot be readily 
submitted to experiment ; but the probability that f{v) should satisfy 
the solution for two dimensions, and fail in that for three, is so slight 
that it may be neglected. 

Our first experiments were with a small iron plate covered with a 
mixture of wax, rosin, and paraffine. By this means, one curve was 
obtained at each heating ; viz., that separating the part of the mixture 
which had melted from that which remained solid. This method can 
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be used in the open air, but is impracticable when the waxed plate is 
covered by a non-conducting material. 

After some further rude work, we constructed the table of which a 
diagram is given below : — 





Two Bunsen burners are enclosed in an iron case, which is itself 
enclosed in a wooden case, surrounding it at a distance of 5 to 7 cm. 
Through the top of the iron case is put a bent copper rod, heated at 
either elbow by one of the Bunsen flames. The wooden case rests on 
a floor fastened beneath the table. The ends of the rod, rising from 
the iron case, and surrounded and held firmly in place by tin cylinders 
closed at both ends, and themselves secured to enclosing wooden pipes, 
project some 5 cm. above the top of the table. 

Around this top is a guard, within which are wooden cleats, by 
means of which to level the non-conducting material. This non-con- 
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ductor is the almost pure silicious earth, dug in Keeiie, N. H., and 
known as "infusorial earth." With it, we were abundautly supplied 
by the kindness of Messrs. J. A. Wright and James H. Wilson of 
Keene. It is in every way suited to work of this kind, being clean, 
free from any appreciable amount of moisture, and an almost perfect 
non-conductor of heat. 

Upon the ends of the copper rod, covered above and below to with- 
in 2 or 3 cm. of the edges, rests a sheet of No. 11 boiler-plate iron, 
1.5 meters long and .9 meters wide. 

The first thing necessary to the success of our experiments with the 
plate was that the head of gas should be constant. We found the 
variation to be insensible, although a change in pressure of a fraction 
of one mm. of water could have been detected. 

The temperature of the air must also be constant. To secure this, 
we were obliged to use the precautions mentioned above. The air for 
the lamps entered at the bottom of the cases, and passed out from both 
by an iron pipe placed within a wooden one. The tin cylinders were 
filled with " infusorial." The variation in the temperature of the air 
was less than one degree during the day. It would have been impos- 
sible to keep the temperature exact to the i^^^j^ of a degree, as is 
reported of Biot. 

No work could be done until the plate reached a final state. For 
this an average of five hours was required. After this time the tem- 
peratures were sensibly constant. 

The large size of the plate was necessary, because observations 
taken near the edge ought not to be relied on, and because the varia- 
tion in temperature was quite rapid, whilst even a slight error in deter- 
mining a point on the curve produced a great change in its equation. 
For example, in one case, in a curve whose loop was about 210 mm. 
in diameter, a change of four millimeters in the position of a given 
point changed m in the equation r^r^"^ = c from — .5 to -[- '^j or the 
equation from 

-— = c to r.rJt = c. 

In another case, a change of three millimeters altered the equation from 

r^r^^ =zc to r^rA = c 

The measurements were made with the galvanometer and thermo- 
pile. Our galvanometer has four coils, whose resistance is 11. ohms, 
and a mirror of 1 meter focal length. Each scale-division is f mm. 
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In order that the electromotive force of the current, and therefore 
the deflections of the galvanometer, should be directly proportional to 
the excess of temperature of the point touched over that of the air, 
a pair of metals must be taken whose neutral point is at a very high 
or very low temperature, and whose relation to each other varies the 
least possible with the temperature. Iron and German-silver satisfy 
these conditions remarkably well. Their neutral point is at 1354*^, 
as deduced from Prof. Tait's table. Some rude experiments of ours 
on the relation between the temperature and the deflections gave quite 
ragged curves ; but they were evidently in every case straight lines. 
The three lines obtained by heating and cooling water were nearly 
parallel. The mean of these experiments showed that one degree of 
temperature corresponds to 4.3 scale-divisions. One microvolt corre- 
sponds to 2.6 scale-divisions, or one scale-division to .39 microvolts. 
We thus find that one degree of change in temperature produces an 
electromotive force of 18.3 microvolts. The value as calculated by 
Prof. Tait's formula is 21.8. 

Our first experiments were with two thermopiles. The method of 
work was this. Each thermopile consisted of one joint, which was 
fitted to a piston moving vertically 3 cm. in a brass cylinder. The 
piston was depressed by a weight resting on a lever, whose fulcrum 
was on the frame around the table, and its end at the thermopile. 
The average pressure applied was 9 kilos. This was opposed by the 
reaction of springs-, which raised the piston when the weight was re- 
moved. The cylinders were placed on carriages, running on wooden 
strips, which crossed the table transversely just above the infusorial. 
These strips and the sides of the table were graduated, and thus we 
got the. (x, y) co-ordinates of the points touched. 

One thermopile was kept in a fixed position, the other was moved 
about to various points where equal deflections were obtained. All 
these last points lie on an isothermal. Great care is necessary that 
the contact with the plate be good, and that the piston move in an 
exactly vertical direction. The results thus found were unsatisfactory, 
because, as before pointed out, an inappreciable error in observation 
produces a very perceptible change in the equation of the curve. An- 
other objection is that the method requires a great expenditure of 
time. 

A second method is more satisfactory. One thermopile was moved 
lengthwise along the plate, and the deflections laid off in a curve. 
The same may also be done transversely, but the curves obtained are 
of little value. 
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The results of our work on " longitudinal curves " are given below. 
Evidently, in crossing the plate, we shall find, in general, four points 
having the same temperature, two about either pole. The isothermal 
points lying about one pole in each of two experiments, obtained by 
averaging several readings of the curves mapped out, are as follows : — 





Isothermal Points. 






Exp. 1. 


Exp. 2. 


Mean. 


488 —658.1 


488 


— 661.9 


488 


— 660. 


495 —644.4 


495 


— 643.4 


495 


— 643.9 


500 _ 632.6 


500 


— 631.1 


500 


— 631.85 


502^ — 626.8 


502i 


— 625.9 


502i 


— 626.35 


505 —622.5 


505 


— 620.6 


505 


— 621.55 


507^-617.4 


507i 


— 615.9 


507i 


— 616.65 


510 —612.9 


510 


— 611.8 


510 


— 612.35 


515 —604.4 


515 


— 602.3 


515 


— 603.85 


517 —600.9 


517 


— 599.4 


517 


— 600.15 


520 —596.9 


520 


— 595.6 


520 


— 596.25 


525 —591.3 


525 


— 588.8 


525 


— 590.05 


527 —588.4 


527 


— 586.6 


527 


■ — 587.5 


530 —584.6 


530 


— 582.5 


530 


— 583.55 


535 —578.7 


535 


— 576.2 


535 


— 577.45 



555 — 555 555 — 555 555 — 555 

The isothermal curves corresponding to these pairs of points are in 
general form like the family to which the common lemniscate belongs. 
Several other curves have much the same form, and therefore it be- 
comes necessary to show that our experimental curve is none of these. 
The following conditions must be satisfied : — 

(1) For a given positive value of r^ there must be only one real 
and positive value of r^, r^ and rg being radii from the fixed points. 

(2) The function must be homogeneous ; else, for the same physi- 
cal problem we shall get as many different solutions as we employ 
units of measure. 

Hence we must exclude trigonometric functions, by reason of their 
periodicity, and, since exponential and logarithmic forms may be devel- 
oped in algebraic series, need consider only these last. 

In this system of coordinates, negative, as well as truly imaginary 
values of r^ and r^, are to be regarded as imaginary. (1) requires us to 



OP ARTS AND SCIENCES. 



135 



exclude any form in which both positive and negative terms enter, as 
also forms where the exponents of the powers of r are both positive 
and negative. 

Three simple equations satisfy these conditions : — 

(1.) r,r, = c. (2.) J- + J- = c. (3.) r^ + r, = c. 

(3) represents the ellipse; and all forms r^ -^ r^ z=^ c are of 
equally little use to us. 

The isothermal points for 1 = c are approximately : — 

488 —627 515 — 595^ 

495 — 620 517 — 593 

500 —615^ 520 — 590 

502^ — 612 525 — 585J. 

505 — 608^ 527 — 583| 

507 J — 605 530 — 580 

510 — 602i 535 — 575 

The isothermal curves for — - ~| r = ^ are much nearer circles, and 

the isothermal points more nearly equidistant from the pole. This 
approximation to a series of circles increases with the degree of the 
equation. There remains to be considered only 2^ r^-^* r^ -*=: c. 
The physical conditions require that the two loops of the isothermal 

curves shall be approximately equal. But, whenever the ratio 

m — X 

in the equation r^"* "*" * rg** ' * = c differs much from 1, the loops are very 

unequal ; and any term in which a ratio of this sort is found introduces 

into the equation 2J r^*^ + * rg** ■" * = c an inequality unbalanced by the 

other terms. 

Assuming now r^r^^ = c, as the equation of our experimental curve, 

the values of p, corresponding to the isothermal points in the " Mean " 

column, are as follows : — 



(1) 


1.097 


(8) 


.988 ' 




(2) 


1.125 


(9) 


.912 




(3) 


1.093 


(10) 


.950 




(4) 


1.072 


(11) 


1.056 


average 1.059 


(5) 


1.065 


(12) 


1.089 




(6) 


1.041 


(13) 


1.098 




(7) 


1.035 


(14) 


1.207 
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Compare now the experimental values before given with the four 
series of values below. 

The third and fourth fail to satisfy condition (2), as this condition 
was discovered too late to change our figures, but they will serve 
equally well for illustration. 

Values corre- in in in in 

668. 

649. 

637. 

630.5 

625. 

620. 

615. 

606.5 

602.5 

598. 

591.5 

588. 

583. 

577. 

It will appear at once that while the first two curves lie very close 
to the experimental curve, and parallel with it, even such close approx- 
imations to the one-term equation as those above cut at an angle. The 
use of residual curves will make the fact still clearer. Equations of 
three, four, or more terms, will yield like results. 

Even if an equation could be found of this kind which was equally 
good with the simple form, it would still be improbable that the sim- 
pler should not be the form preferred by nature. 

Since then the temperature is shown experimentally to be constant 
along the family of curves r-^r^ =z c, a function f(v) can be found, 
the same for all bodies, whose derivative in any direction in a given 
body, when multiplied by a constant quantity depending on the nature 
of the body, measures the flux of heat in that direction. 

* Computed. t Approximate. 



sponding to 


r^r^l'^ = c* 


^^^^1.075 == c* 


rir^ + rir^'- 


488 


657.1 


658.5 


662. 


495 


640.8 


641.8 


613. 


500 


630.6 


631.3 


632. 


502^ 


625.8 


626.4 


625. 


505 


621.2 


621.8 


620.5 


507i 


616.8 


617.3 


616.5 


510 


612.6 


613. 


612. 


515 


604.6 


604.8 


602. 


517 


601.5 


601.8 


600. 


520 


597.1 


597.3 


595. 


525 


590. 


590.2 


588. 


527 


587.3 


587.4 


585.5 


530 


583.4 


583.5 


581. 


535 


577.1 


577.15 


575.5 
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The solution for one dimension ' ^ = , as well as that for 

two, admits of experimental verification. 

Let a long rod, covered with non-conducting material, except at the 
ends, be heated at the points A and B. Then for A : f(i\) = k^r^^ -|- 
by For B: f{v^ z=z k^r^ + ^2» ^i ^"^ ^2 being the distances from A 
and B measured on the rod. 

f{v) =f(v,) +/(r,) = Vi + V2 + ^ + ^2 = Vi + V2 + ^• 

Let k^ = Ic^ fvz=: k^(r^ + r^) + L 

f{v) between A and B is constant, and may be represented by a 
straight horizontal line. Whenever /(i;) has a constant value, v has a 
constant value, and conversely. Therefore, v, as well as /(?;), must be 
represented between A and ^ by a straight line. If A and B are not 
heated exactly alike, the line will still be straight, but inclined to the 
horizontal. 

We have not yet tried the experiment, but intend to do so, that 
there may be the greatest possible assurance in regard to the existence 
of f(v), seeing that it exists in every case which can be experimentally 
tested. 

The second portion of our work is the determination oi f{v). For 

a rod heated to the final state, it must satisfy the condition ''• ^ = 0, 
ovf{v) = Ax -\- B. 

For a plate: ^^^ + ^-^^ = 0, ov f{v) = G+D log r (chang- 
ing to polar co-ordinates). 

For a solid: -^ + -^^ + — _ = O,or/(r) = ^+ -. 

This function cannot be v itself. For, in the case of a rod, heated at 
one point, let/v = v. 

V = Ar-\-B. 

When r = l v z=z A -\- B ^(X) 

„ r = V = j5=:[=Q0 

„ 7' =z CO v = Aoo-\-B=Oj which is impossible. 

For a plate 

f{v) =v= (7+ Dlogr 

W^hen r=l i; = ^=|=oo 

„ r = V z=. G -\- ( — 00 ), which again 

18 contrary to the physical conditions. 



since by physical conditions 
the temperature must here 
be finite. 
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For a solid 



When r=cc v == .\ E = 



F 

V = —- 
r 



F 

„ r = V =z — ^=z Qc , which again is impossible. 

The same thing may be proved as follows. Assume, for the sake of 
argument, that the flux is — c—- for the direction x, — c-j- for y, 

— c— for z. Then for all homotjeneous bodies. 



dv _ /d^ \ "Hi \ iiax 

dt~ Xdx'^ "^ df "^ dz^y 



This is the only condition that a function v must satisfy, in order to 
represent an actual possible case. Any function that is a particular 
solution of this equation may represent an actual distribution of heat. 
First, liv =z cp(x, y,z, t) represent the temperature throughout a 
body, it cannot have any true maxima and minima for x, y, z. This 

is evident, since the conditions of a maximum or minimum are -4- = 0, 

-4- = 0, -J- = 0, which by hypothesis cannot be zero without mak- 
ing the flux = 0. No point can be hotter or colder than the points 
around it, and there not be a flux to or from the point. The physical 
conditions forbid the mathematical condition of the existence of maxi- 
ma and minima. There must be points hotter than points around 
them, and therefore they must be shooting points, and not maxima or 
minima. 

Secondly. If t;^ and v,^ are two particular solutions of the Partial 
Differential Equation, their sum is also a solution, and therefore cor- 
responds to an actual distribution of heat, in which the temperature of 
any point is equal to the sum of the temperatures, which it would have 
under the conditions represented by v^ and v^^ It will now be easy to 
show 

Thirdly. The points of hottest temperature, when the solution is v, 
fall in exactly the same places in the body as the hottest points of the 
two solutions v^ and v^ taken jointly. That is to say : if there are w^ 
points of highest temperature when the solution is v^, and 71^ when the 
solution is v^i there will be n^ -\- n^ points of highest temperature when 
the solution is v^ -^^ v^ = v, unless v^ and v^ have some hot points in 
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common, or else a mkiimum of v^ or v^ corresponds with a maximum 
of Vg or v^. Here special investigation is necessary for each particu- 
lar case. 

Observe carefully that, if the functions v, r^ and v^ attained true 
maxima or minima at the hot points, instead of coming to a point, the 
proposition would not be true. 

That the proposition is true, in regard to shooting points, can be seen 

thus : ^ = tan t ; tan r changes sign instantaneously at a shooting 

point, and hence, for a small (infinitely small) change in a;[6?x], ^ 

has a finite change of value, and d^-j- = finite, dx =. ^. There- 
fore 

dxV finite 

^-m = 1 =00 



dx 



This is the condition for a shooting point. If, at the same point, -4- 

d^\3 
as well as -— is infinite, there is a cusp. If also v = x , there is 

an asymptote. 

At the points of highest temperature -~ = oo , ~^ = oo , — 

= oo . 

Now — --- = ~— H 

dx^ dx'^ ~ dx'^ 

and any values of ic, y, z, which make either term in the second member 
infinite, will make the whole member infinite, and insure a shooting 
point for v as far as the direction x is concerned, wherever there is one 

for either v^ or v^. The like is true of y and z. If -^-^ and -^ - become 

dx^ dx'^ 

00 with opposite signs, there will be an ambiguity which can be got 

rid of by the determination of the indeterminate quantity (oo — oo ). 

Using the results of both (2) and (3), it is easy to arrive at an 
absurdity in almost any case considered. For instance, heat to a con- 
stant temperature all the points upon a circle, marked out upon a large 
metal plate protected from surface radiation. Then not only will the 
point at the centre of the circle soon become hotter than any point upon 
the rim, in which case the heat must flow from within out, but it must 
finally be infinitely hot, which is absurd. 

Our experiments have not proceeded so far that we can determine 
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f(v) from them. That determination must form the subject of a sub- 
sequent paper. We will simply present here a specimen of our exper- 
imental results for both the rod and the plate, heated at one point 
under the conditions above specified ; the values in column " v " being 
galvanometer deflections : — 

Eod. Plate. 



r 


V 





342 


100 


257 


200 


196 


300 


150 


400 


111 


500 


87 


600 


67 


700 


54 


800 


43 


900 


33 


1000 


29 



r 



V 

178 


50 


108 


100 


83 


150 


61i 


200 


53 


300 


37 


400 


281- 


500 


20i 


600 


15 



From these series we are to determine v as a function of r. If v 
= qj{r). Then for a rod, fcp{r) =z Ar -\- B; for a plate, yg)(r) = 
0+i>log r. 

Permit us here to express our great gratitude to your Academy for 
the generosity and liberality which have supplied us with the means for 
carrying on our experiments. 

Harvard University, Oct. 10, 1877. 



